Orbitally driven spin-singlet dimerization in S=l La4Ru20io 
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Using x-ray absorption spectroscopy at the Ru-Z/2,3 edge we reveal that the Ru*"'' ions remain in 
the S=l spin state across the rare 4d-orbital ordering transition and spin-gap formation. We find 
using local spin density approximation -f Hubbard U (LSDA+U) band structure calculations that 
the crystal fields in the low temperature phase are not strong enough to stabilize the S=Q state. 
Instead, we identify a distinct orbital ordering with a significant anisotropy of the antiferromagnetic 
exchange couplings. We conclude that La4Ru20io appears to be a novel material in which the 
orbital physics drives the formation of spin-singlet dimers in a quasi 2-dimensional 5=1 system. 

PACS numbers: 71.20.-b, 75.25.+Z, 75.30.Et, 78. 70. Dm 



One of the most intriguing aspects of transition metal 
materials is the wide variety and richness of their physical 
properties Q. Although conceptually clean and beauti- 
ful, theoretical simplifications in terms of a Hcisenbcrg 
model or a single band Hubbard model turn out to be in- 
adequate |3| ■ It now becomes more and more clear that a 
full identification of the relevant orbital and spin degrees 
of freedom of the ions involved is needed to understand, 
for instance, the colossal magneto-resistance behavior in 
the manganates 0, 0, 111 j magnetization reversals and 
metal-insulator transitions in early transition metal ox- 
ides H S, 0, E3 , as well as the formation of spin-gaps 



in non-l-dimensional S=\ systems |llUl2l ll3lll 

Very recently Khalifah et al. synthesized the semi- 
conducting quasi-2-dimensional La4Ru2 0io compound 
and discovered that this system undergoes a strong first- 
order structural transition at Ts=160 K (see Fig. 1), 
accompanied by a rare 4d-orbital ordering and spin-gap 
opening. Their interpretation of these phenomena was 
that the Ru'*+ ion transforms from the usual ^231^231 
low-spin state with 5=1 to an t^g^t^g^ 'ultra-low' spin 
state with S=0, caused by a sufficiently strong crystal- 
field splitting (CFS) of the Ru Ad-t2g levels due to the 
lattice distortion below Tg. 

However, already soon after that, it was also hypoth- 
esized by Khalifah et al. based on unpublished standard 
band structure calculations that a chemical bond may 
be formed associated with the orbital ordering. Here we 
report on an x-ray absorption spectroscopy (XAS) study 
in which we reveal that the Ru^+ ions remain in the iS'=l 
spin state across Tg. This directly points to the possi- 
bility that La4Ru20io is in fact a novel system in which 
the spin-gap opening is due to a singlet dimer formation 




FIG. 1: (Colour online) The main block of the low- 
temperature (20 K) triclinic crystal structure of La4Ru20io 
shown in the crystallographic ab plane: The 2-dimensional 
Ru-O network consisting of distorted RuOe octahedra extends 
along the c-axis (not shown, pointing into paper plane) and b- 
axis with corrugation. The Ru-O bond-lengths are marked in 
unit of A, and those out of the aft-plane (not shown) are 2.046 
and 2.057 A for Rul, and 2.046 and 2.082 A for Ru2. The 
spin-singlet dimers are marked by black solid bars. The local 
orthogonal coordinate system (xyz) is used in our band cal- 
culations with z parallel to c and y (x) along the short (long) 
Rul-Ru2 direction. The high-temperature monoclinic struc- 
ture (not shown) above rs=160 K has equal Ru-Ru distances 
along the x and y directions: Rul and Ru2 are equivalent. 



in a non-l-dimcnsional and S>^ material. We find using 
LSDA+U band structure calculations that the distinct 
orbital ordering involves a significant anisotropy of the 
antiferromagnetic exchange couplings, indicating indeed 
the formation of Ru(S'=l)-Ru(S'=l) singlet dimers. 

Floating zone crystals were grown in an NEC SC- 
M15HD image furnace using rods with a 1:1 or 4:5 ra- 
tio of La:Ru which had been pre-reacted and sintered in 
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air at 1250°C. A small (1-2 atm) overpressure of oxygen 
aided the growth, and the power was dynamically in- 
creased during the run to compensate for absorption by 
the copious amounts of evaporated Ru. Sizeable crystals 
could only be obtained using seed crystals. X-ray diffrac- 
tion confirmed both the macroscopic phase purity and 
the universal presence of two twin domains. The XAS 
measurements were performed at the Taiwan NSRRC 
15B beamline, equipped with a double Si(lll) crystal 
monochromator delivering photons from 2 keV and up. 
The spectra were recorded using the total electron yield 
method in a chamber with a base pressure in the low 
10"^*^ mbar range. Clean sample areas were obtained by 
cleaving the crystals in-situ. The photon energy resolu- 
tion at the Ru 2^2,3 edges (/ii/w2.9 keV) was set at 0.6 
eV. Strong polarization dependent 0-K XAS spectra 0| 
verify the high sample quality. 
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FIG. 2; (Colour online) . (a) Experimental RU-L2, 3 XAS spec- 
tra of La4Ru20io measured at 220 and 110 K, i.e. above and 
below Tg — IGQ K. (b) Theoretical simulations for the Ru*"^ 
ion in the S=l high-temperature phase at 220 K, S=l low- 
temperature phase at 110 K, and 5=0. The inset shows O-K 
spectra measured at 200 and 85 K from fl^ . 

The top curves of Fig. 2 depict the Ru-L2,3 XAS spec- 
tra of La4Ru20io taken at 220 K (black line) and 110 
K (red line). The spectral line shapes depend strongly 
on the multiplet structure given by the Ru Ad-Ad and 
2p-Ad Coulomb and exchange interactions, as well as by 
the local CFS and the hybridization with the O 2p lig- 
ands. Unique to XAS is that the dipole selection rules 
are very effective in determining which of the 2p^4d"+^ 
final states can be reached and with what intensity, start- 
ing from a particular 2p^4c?" initial state (n=4 for Ru^+) 
|l7l Il8l | . This makes the technique extremely sensitive to 
the quantum numbers of the initial state [l9|, |20j . 

The essence of the Ru-L2,3 XAS spectra in Fig. 2 is 
that there is only a very small change across T^, suggest- 
ing that the local electronic and spin state of the Ru^+ 
ion in the high temperature (HT) and low temperature 
(LT) phases are quite similar. This seems surprising in 



view of the fact that we observed considerable modifica- 
tions in the O-K XAS spectra in going from 200 K (black 
line) to 85 K (red line) as shown in the inset of Fig. 2 
[lij . Since these O-K spectra are known to be sensitive 
to band structure effects ^3 1 their modifications are fully 
consistent with the strong changes in the crystal struc- 
ture as seen in neutron diffraction [isj . confirming once 
again the good quality of our samples. 

To extract quantitative information on the CFS and 
spin state from the Ru spectra, we have performed simu- 
lations using the successful configuration interaction clus- 
ter model |l7lll8ll21| . The calculations have been carried 
out for a RuOe cluster in the proper HT and LT local co- 
ordination using the XTLS 8.3 program Parameters 
for the multipole part of the Coulomb interactions were 
set standardly at 80% of the Hartree-Fock values [2l| . 
while the monopole parts {Uddi Upd) were taken from 
Ca2Ru04 [2II23. The O 2p - Ru Ad charge transfer en- 
ergy was estimated from LDA calculations (sec below), 
and the O 2p - Ru Ad transfer integrals and their distance 
dependence from Harrison's relations The local CFS 
parameters are to be determined from the comparison 
between the simulations and the experimental spectra. 

The bottom curves of Fig. 2 show the simulations for 
both the HT and LT phases. We found optimal fits (black 
and red lines) if the Ru-4d xz orbital is set at about 100- 
150 meV (HT) and 200-300 meV (LT) lower in energy 
than the essentially degenerate (within 50 meV) yz and 
xy orbitals jl^. These numbers refer to total energies 
calculated including the CFS and covalcncy but without 
spin-orbit interaction. Important is that the cluster cal- 
culation indeed finds the iS'=l state for both the HT and 
LT phases, which is a direct consequence of the fact that 
in both phases the xz orbital is essentially doubly occu- 
pied while the yz and xy are each singly occupied. We 
also have carried out calculations for the Ru ion in the 
artificial S—0 state by changing the CFS parameter such 
that the xz orbital lies above the degenerate yz and xy 
orbitals [2^. As shown in Fig. 2, the simulated spectrum 
(green curve) disagrees with the experiment. So we can 
rule out that the spin-gap opening is due to a local spin 
state transition [15|. 

To confirm our XAS-derived conclusions and, more im- 
portantly, to get in-depth understanding of the nature of 
the spin-gap state below T,, we performed systematic 
LDA and LSDA-I-U band structure calculations jl^, by 
using the full-potential augmented plane waves plus lo- 
cal orbital method [l^. We took the neutron crystal 
structure data at 20 K and 298 K The muffin-tin 

sphere radii are chosen to be 2.8, 2.0 and 1.5 Bohr for La, 
Ru and O atoms, respectively. The cut-off energy of 13 
Ryd is used for plane wave expansion of interstitial wave 
functions, and 120 k-mesh for integrations over the Bril- 
louin zone. U=3 eV and Hund's rule exchange Jij=0.5 
cV (Ue//=2.5 eV) are used for Ru Ad electrons, which 
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FIG. 3: (Colour online). LDA density of states of the Ru 
4d orbitals in the high-temperature (HT, upper panel) and 
low-temperature (LT, middle and lower panels) phases. The 
insets show a close up of the t2g {xy, xz, and yz) states in the 
vicinity of the Fermi level set at eV. 

Fig. 3 shows the Ru Ad density of states (DOS) cal- 
culated using the LDA for the non-magnetic (NM) state. 
The inset shows a close-up of the t2g levels and a calcu- 
lation of the first moments supports the XAS analysis: 
for both the HT and LT phases, the xz orbital lies low- 
est and the splitting A,^^ between the higher lying yz 
and xy orbitals is less than 50 meV. It is the magnitude 
of this App. relative to in which determines the spin 
state of the Ru*+ ion. To make a crude estimate: the 
S=\ state {xz^^yz^ xy^) carries the Hund's stabilization 
energy of 3J//, whereas the 5=0 state {xz^^yz"^^) has a 
total stabilization energy of 2Jh plus A^^^. Assuming 
that 0.5 eV is a reasonable estimate for Jjy, one must 
expect that a A^^ of 0.05 eV is far from sufficient to 
obtain the 5=0 state. Obviously this is what the XAS 
experiments have revealed. Moreover, our LDA calcula- 
tions also find that the ferromagnetic (FM) as well as the 
antiferromagnetic (AF) solution are more stable than the 
NM one, giving further support that the S=0 is unfavor- 
able. This was also recently found by GGA calculations 
of Eyert et al^. 

The LDA results as shown in Fig. 3 predict a metallic 
state for both the HT and LT phases, and this is in strong 
disagreement with the observed semiconducting behavior 
with a gap of about 0.3 eV One should note how- 
ever, that the Fermi level is located in a narrow Ru t2g 
band with no more than 1.5 eV width. This signals that 
modest electron correlation effects at the Ru sites will 
already be able to turn this material into a Mott insu- 
lator. We therefore set out to do LSDA-I-U calculations 
for the LT phase and found that a band gap of about 0.5 
eV is indeed opened in the Ru t2g band as can be seen 



from Fig. 4a. Within the LSDA-I-U mean-field approach, 
the lowest state of this insulator is AF and is labeled as 
AFxyz in Table I to indicate the AF alignment with the 
nearest neighbors along the a;, y, and z directions (Fig. 
1). The Ru*^ spin moment inside the mufhn-tin sphere 
is about 1.2 /is, confirming the XAS result that the spin 
state is S=\ but not 5=0. 

It is important to look now at the orbital character of 
the AYxyz solution. Fig. 4b and 4c show the orbitally 
resolved DOS of the AF aligned Rul and Ru2 ions. One 
can clearly see for each of the ions that the xy and yz or- 
bitals (with the spins parallel) are singly occupied while 
the xz are doubly occupied. This double occupation is 
due to the Ru-0 bonds being elongated along the x di- 
rection, see Fig. 1. We thus find an orbital ordered state 
which is different from the originally proposed doubly 
occupied xy and xz (5=0) state [l5| . 

As a result, the half- filled xy and yz orbitals are mag- 
netically active. To explain the formation of the spin 
gap in the LT phase, it is crucial to identify the relevant 
exchange interactions in this system. We therefore have 
calculated other magnetic configurations: the FM and 
two more types of AF solutions, namely AF^^j, and AF^^. 
The 2xLa4Ru2 0io unit cell is used for all solutions, ex- 
cept for AFj,2 where the 4xLa4Ru20io unit cell is taken 
with the doubling along the b direction (sec Fig. 1). We 
also studied the NM solution, and confirmed that this is 
much higher in energy than the AF^yz^ by 775 meV per 
formula unit (Table I). The NM solution is metallic, the 
FM half-metallic and all AF insulating. The Ru*+ spin 
moment is 1.2 ± 0.1 ^lb for all magnetic solutions. 

The relative energies of the different magnetic states 
allow us to estimate the exchange constants Jx, Jy, and 
Jz along the x, y, and z directions (Fig. 1), respectively. 
With the FM being 130-150 meV higher in energy than 
the AF solutions, we thus have very large AF exchange 
interactions in this system. As listed in Table I, we find 
Jj,=1.5 meV, Jj,=65.5 meV and Jz=4.5 meV. The sig- 
nificant anisotropy is related to the fact that both the xy 
and yz orbitals contribute to the exchange coupling along 
the y direction having the short Ru-O-Ru distance, while 
only the xy or yz orbital contributes to the x or z di- 
rection, respectively, having the long Ru-O-Ru distance. 
Hence we can consider the LT phase of La4Ru20io as 
practically consisting of strongly coupled Ru-Ru dimers 
with weak intcr-dimer coupling, or at most as two-leg lad- 
ders with Jrung~Jy=G5.5 mcV and Jieg=Jz=A.b meV, 
with weak interladdcr coupling of 1.5 meV. Those rather 
well isolated Ru-Ru dimers or rungs will have the singlet 
ground state [sij . This explains naturally the appear- 
ance of a spin gap in the LT phase, which, according to 
our calculations should be about 60 meV, in reasonable 
agreement with the experimentally observed value of 40 
meV |l5j . The dimer character of the spin gap seems 
also to agree with the results of the single-crystal ncu- 
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FIG. 4; Density of states (DOS) of the La4Ru20io in the antiferromagnetic low temperature phase calculated using LSDA+U. 
(a) Total DOS per formula unit, the 4d states of the two inequivalent Ru'*^ ions, and the 2p states of all the ten oxygens; (b) 
orbitally-resolved 4d states of Rul and (c) of Ru2. The solid (dashed) curves denote the up (down)-spin channels. 



tron scattering [32 



TABLE I: LSDA+U results for the total energy (per formula 
unit and relative to the lowest solution) and band gap in the 
low-temperature phase of La4Ru20io, calculated for the non- 
magnetic (NM), ferromagnetic (FM) and three types of an- 
tiferromagnetic (AF) solutions. The exchange constants are 
found to be J3;=1.5 meV, J^=65.5 meV, and Jz =4.5 meV. 



state 


exchange 


energy (meV) 


gap (eV) 


NM 




775 




FM 
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152 




AF 

-^J- xyz 


Jx Jy 2Jz 





0.5 


AF 
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To summarize, XAS measurements revealed that the 
Ru^+ ions in La4Ru2 0io remain in the S=l spin state 
across the structural phase transition and spin-gap for- 
mation. LSDA-f U calculations provided support for this 
finding and identified the distinct orbital ordering accom- 
panying the structural transition. Crucial is that with 
LSDA-fU we were able to estimate the inter-site anti- 
ferromagnetic exchange interactions and found them to 
be highly anisotropic. This brought us to the conclu- 
sion that the spin-gap opening is due to the formation 
of Ru-Ru singlet dimers. Such a transition is rather un- 
usual since La4Ru20io is a 2-dimensional S=\ system; it 
is largely driven by orbital ordering which amplifies the 
importance of orbital physics in correlated systems. 
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